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Abstract — A Rotating Modulator (RM) is one of a class of 
techniques for indirect imaging of an object scene by modulation 
and detection of incident photons. Comparison of the RM to 
more common imaging techniques, the Rotating Modulation Col- 
limator and the coded aperture, reveals trade-offs in instrument 
weight and complexity, sensitivity, angular resolution, and image 
fidelity. In the case of a high-energy (hundreds of keV to MeV), 
wide fleld-of-view, satellite or balloon-borne astrophysical survey 
mission, the RM is shown to be an attractive option when coupled 
with a reconstruction algorithm that can simultaneously achieve 
super-resolution and suppress fluctuations arising from statistical 
noise. We describe the Noise-Compensating Algebraic Recon- 
struction (NCAR) algorithm, which is shown to perform better 
than traditional deconvolution techniques for most object scene 
distributions. Results from Monte Carlo simulations demonstrate 
that NCAR achieves super-resolution, can resolve multiple point 
sources and complex distributions, and manifests noise as fuzzy 
sidelobes about the true source location, rather than spurious 
peaks elsewhere in the image as seen with other techniques. 

Index Terms — Gamma-ray Astronomy, Image Reconstruction, 
Multiplexing, Scintillation Detectors. 

I. Introduction 

IMAGING hard x-ray and gamma-ray photons (hundreds 
of keV to MeV) cannot be accomplished using focusing 
techniques, as with lower energy photons. A concise overview 
of several techniques is given in HI. In one class of methods, 
incident photons are spatially or temporally modulated before 
detection. The recorded data are not a direct representation 
of the object scene, and so additional steps are required to 
deconvolve this information with a pre-determined instrument 
response function (i.e., system matrix). 

A variety of deconvolution techniques have been developed 
across a wide range of applications, each employed for its 
demonstration of computational speed, noise suppression, re- 
solving power, or fidelity; classes include statistical, algebraic, 
and "ad-hoc" algorithms. Algebraic techniques attempt to 
solve directly for the unknown image. Since this class requires 
convergence of the reconstruction to the data, noise in the data 
can be amplified and cause spurious peaks to arise. As we will 
show, however, an algebraic technique has a distinct advantage 
(the ability to achieve "super-resolution") that makes it entirely 
suitable for the rotating modulator, and with an appropriate 
non-linear step, statistical noise may be adequately suppressed. 

II. Rotating Modulation 

A Rotating Modulator (RM) |2), is one of a class of 
temporal modulation imagers. It consists of a single grid of 
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opaque slats spaced apart by equally-wide slits, suspended 
above a small array of circular non-imaging detectors (Fig. 
[]J). The detector diameters are equal to the slat/slit widths. 
The grid rotates, periodically blocking the transmission of 
incident photons from the object scene onto the detection 
plane. For each detector, this time-dependent shadow generates 
a characteristic time history of counts for the entire length 
of the exposure-varying from to 100% transmission-which 
is then folded modulo the rotational period of the grid. The 
profiles are combined with a collection of pre-calculated 
profiles to generate a cross-correlation image. Deconvolution 
of the cross-correlation with the instrument response produces 
the final image reconstruction. 

To best understand the characteristics of the single-grid 
RM in the scope of the modulation class of imagers, it is 
useful to consider the better-known temporal imager, the two- 
grid Rotating Modulation Collimator (RMC) £Q, 0, and the 
traditional spatial modulator for high-resolution gamma-ray 
imaging, the coded aperture J6), Q. The latter incorporates 
a mask of opaque and transparent pixels (typically an equal 
number of each, resulting in a 50% mask transmission), which 
spatially modulates incident photons onto a position-sensitive 
detection plane beneath (Fig. [TJ)). Angular resolution, 58, is 
defined by the ratio of mask pixel size, a, to mask/detector 
separation, L: 58 = a/L. The detection plane should have 
spatial resolution half the mask pixel size or smaller (to satisfy 
the Nyquist condition) in order to resolve a shadow well 
enough to obtain the intrinsic resolution of the instrument. For 
a given detector pixel size, the maximum achievable resolution 
is limited due to a finite sampling of the shadow pattern. 
The resolving power may be increased only by further sub- 
dividing the detection plane. For a high-sensitivity, large-area 
instrument, this results in a high number of readout channels, 
increasing cost and complexity. (In the case of the CASTER 
and EXIST 0, (lO) designs for a satellite-based Black 
Hole Finder Probe gamma-ray telescope, the requirement for 
10 minute-of-arc-scale angular resolution and meter-squared- 
scale sensitive area leads to in excess of 10 4 - 10 7 readout 
channels.) 

An RMC (Fig. [lj) consists of two offset concentric grids 
of equal-width slats and slits. One grid sits just above a single 
non-imaging detector, while the other grid is a distance L 
overhead. As the two grids rotate in tandem, incident photons 
from off-axis sources within the field of view (FOV) will 
be periodically shadowed. A characteristic count profile is 
recorded by the detector, and then used for image recon- 
struction. As with the coded aperture, the geometric angular 
resolution is defined by the ratio of slit width to mask-detector 
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(a) Rotating Modulator 




(b) Coded Aperture 




(c) Rotating Modulation Collimator 

Fig. 1. Basic structure of three multiplexing imagers. All instruments consist 
of a detection plane situated below a suspended mask which spatially or 
temporally modulates incident photons. 



separation, and so grids with small spacing can be constructed 
to provide excellent angular resolution. Due to the continuous 
nature of the measurement vector, the observation profiles 
may be sampled at higher frequency with no increased cost 
or complexity. With an appropriate reconstruction algorithm, 
resolution can be enhanced beyond the geometric resolution, 
although at lower imaging efficiency (see discussion in Sec. 
[V) . At high photon energies, the combination of large grid 
thickness and small slat spacing limits the maximum ob- 
servable source angle due to mechanical collimation. With 
only a single detector, the RMC is a much simpler system 
than a coded aperture, though in practice, a single large- 
area detector would require multiple readout channels. The 
addition of a thick second grid adds to instrument weight and 
reduces the average transmission to 25%. Furthermore, image 
reconstruction generates "mirror" sources, and the RMC is 
insensitive to sources directly overhead. 

The RM may be considered a trade-off between the coded 
aperture and RMC, featuring temporal modulation and modest 
spatial resolution. It features a simple readout system, 50% 
mask transmission, and relatively low weight. In analogy to 
the coded aperture and RMC, the intrinsic angular resolution is 
defined by the ratio of slit width to mask-detector separation. 
Since this width is set equal to the detector diameters in order 
to maximize sensitivity (by increasing the contrast in the count 
profiles), the resolution is effectively limited for a given mask- 
detector separation and detector diameter. Thus, for the RM 
to be a competitive alternative to other common multiplexing 
instruments, a reconstruction technique must be capable of 
resolving images beyond the intrinsic resolution, i.e. it must 
achieve "super-resolution." The RM then compares favorably 
to a coded aperture with respect to its ability to sub-sample the 
observed data and go beyond the instrinsic resolving power 
and favorably to an RMC at high energies due to a wider 
achievable FOV and the weight of only a single mask. 

While the RMC has seen limited use in gamma-ray 
imaging,-e.g. Ariel-V El (3 - 7 keV), SAS-3 [12J (1.5 - 60 
keV), RHESSI £13] (3 keV - 17 MeV)-the authors are aware of 
no RM previously used in astronomy. The WATCH experiment 
JT4ll was a self-described single-grid RMC, although its design 
was perhaps closer to that of an RM with some key differences 
(including slat detectors instead of circular detectors). Other 
than a prototype RM recently constructed in our laboratory 
at Louisiana State University (LSU) [15], we are aware of 
only one other RM in development [16|. Our intention has 
been to demonstrate the feasibility of the RM design as a 
solution to gamma-ray imaging with wide FOV and good 
sensitivity and angular resolution. This requires the ability 
of an image reconstruction algorithm to resolve beyond the 
instrinsic resolution limit of the instrument. We will show that 
a novel reconstruction algorithm based on an algebraic solution 
can produce images with resolution several times better than 
that defined by the instrument geometry. Furthermore, multiple 
and complex source distributions can be well represented, no 
object information is required a priori, and background and 
counting noise can be adequately suppressed. 
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III. The Imaging Problem 

The data observed by a multiplexing instrument, 0(m), may 
be defined by the counts observed in each data bin m: 



Center Detector 



0(m) = P(m, n)S(n) + B(m) + N(m), 



(1) 



where S(n) is the object scene, B(m) is a background offset, 
N(m) is the measured noise from scene and background 
combined, and P(m, n) is the instrument response, a function 
that transforms the scene information from image space, n, to 
data space, m. In practice, a background subtraction removes 
the B(m) term (though the background noise remains), and 
so it is absent from the formulas henceforth. In matrix form, 
Eq. [T] is written more succinctly as 



0„ 



(2) 



Solving this equation for the object scene S n (i.e. decon- 
volving O m and N m with the instrument response P mn ) is 
the goal of image reconstruction. If P mn were non-singular, 
then the object scene could be written 



(3) 



In general, however, P is neither square nor has an inverse. 
Even if P were non-singular, numerical uncertainties make 
this method impossible in practice. For these reasons, a 
precursor to most techniques is to pre-condition the matrix 
P by multiplying both sides of Eq. [2] by the transpose of the 
instrument response: 



(4) 



or more simply, 



H , — P' C _|_ AT 1 
<-^n' — 1 n'n u n J, n ' ) 

where the cross-correlation image is C n > = P^, m O. 

van 



(5) 
n , and 

the point-spread function (PSF) is P' n , n = P^, m P mn . By 
definition, P' is square-symmetric and will tend towards being 
diagonally-dominant for local PSFs that resemble a two- 
dimensional delta function. For this reason, Eq. [5] is a more 
suitable system to solve than Eq. [2] The unknown noise term, 
N n i — P^, m N m , complicates the solution. The deconvolution 
technique employed must therefore compensate and correct for 
noisy data, while reconstructing true sources. 

IV. RM Instrument Response 

As the RM grid rotates, each detector d produces a time 
history of counts, 



O d (t) = J2Pd(t,n)S(n) + N d (t). 



(6) 



The instrument response Pd{t, n) is a collection of pre- 
calculated count profiles for all possible n source locations 
within the object scene. These profiles are determined analyt- 
ically for an RM with slat width a by l3l 



Pd(t,n) 
1 - 7.F I 1 



r(n) 



cos (£(t) + Co] 



mod 2-1 



(7) 




Fig. 2. Deterministic count profiles from three separate detector locations 
for a single point source. The horizontal scale is grid rotation phase, while 
the vertical scale is unit source intensity. 



where £(t) is the angular orientation of the grid at time t. 
The parameter r(n) is the distance in the plane of detection 
between the detector and the axis of the projected grid shadow 
for a source located at scene location n. If the RM detectors are 
stationary, r(n) remains constant with time. For a grid height 
L, this distance is defined in terms of the source azimuth <f> 
and zenith 8 by 



r(n) = \J (.tq + L tan I 



COS ( 



+ (y + Ltanf5sin0) 2 , (8) 



where xq and j/o are the detector position. The instrument 
is assumed to be fixed relative to the sky; however, if the 
instrument moves with a known behavior, a conversion can 
be made to move these coordinates to a global reference 
frame. (Such an operation would be necessary for a balloon- 
borne or satellite instrument.) The function F(r) describes the 
percentage of circular detector area covered by a shadow that 
has moved a fraction of the total diameter, r, across its face: 



F(t) = - cos _1 (l - 2r) - 2 (1 - 2r)(r - 

7T 7T 



(9) 



The coefficient 7 accounts for transmission through the mask 
and is dependent on the mask thickness h, density p, and mass 
absorption coefficient a for a particular photon energy, as well 
as the angle of photon incidence (i.e. zenith): 



7 = 1 — exp 



hpa 
cos 8 



(10) 



It is important to note that Eq. [7] is an approximation only; 
it ignores photons that "clip" the slat edges. (Description of 
an advanced characteristic formula to extend these results by 
including edge clipping effects will be presented separately.) 
Fig.|2]shows count profile examples of the same source viewed 
by detectors at three different positions. 

The transformation is made from the temporal observation 
to image space by a cross-correlation of the data with the 
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Fig. 3. Cross-correlation image surface mesh of a single point source 
centered in the FOV. The concentric "ring" nature of the PSF is apparent. 
(Axes in degrees) 



transpose of the instrument response summed over the detec- 
tors, 



Substituting 0(t) from Eq. [6] into Eq. 11 we see that the 
cross-correlation image is the object scene correlated with a 
single PSF: 

C{n) = ^P'{ri,n)S{ri) + N(n), (12) 



where 



P'(n',n) = Y,J2P d (t,n')P d (t,; 



(13) 



Conceptually, an element of P'(n',n) may be thought of as 
the relative brightness of pixel n' for a point source located at 
location n (or vice-versa, due to symmetry). Eqs. [12] and [13] 
are analogous to the dirty map and dirty beam, respectively, 
in radio astronomy. For the RM, the PSF appears as a central 
peak with concentric "ring" sidelobes (Fig. [3}, which may be 
described mathematically by a zeroth order Bessel function 

urn 

V. RM Sensitivity 

While a detailed examination of RM sensitivity is beyond 
the scope of this paper, it is necessary to present the formula 
for the signal-to-noise ratio (SNR) of a reconstructed image to 
compare the three discussed multiplexing imagers, particularly 
for the case when super-resolution is achieved. For an RM, an 
RMC, or a coded aperture, the statistical significance of an 
observation at location n of the raw image (Eq. [5]l is found by 
dividing the value in that element by the standard deviation in 
that element |2|. For a total source rate S(n) as measured by 
all detectors in the absence of modulation and total background 
rate B, the SNR is approximated for the case where B 3> S(n) 
as 



where T is the source exposure time, and r\ describes the 
increased sub-sampling of the object scene beyond the geo- 
metric resolution 86. Each factor T) of super-resolution will 
therefore result in a decrease in the instrument efficiency 
of rf 2 - According to Eq. 14 the SNR of a multiplexing 



instrument observation is related to a non-imaging instrument 
by the square root of the variance of the instrument response. 
(In units of SlT/B) 1 ^ 2 , the SNR of an observation by a non- 
imaging instrument is 1.) 

For a coded aperture, the mask pattern is composed of 
ones and zeroes, which represent open and closed elements, 
respectively. Coded patterns are specifically chosen with stan- 
dard deviation 0.5 for all n, the maximum deviation possible 
with an instrument that blocks 50% of the incident photons. 
RMC sensitivity can likewise be calculated, but due to lower 
throughput and non-ideal instrument response, the average 
deviation (and thus statistical significance) is approximately 
0.15. Unlike the coded aperture or non-imaging instrument, 
the sensitivity for an RMC is not uniform across the sky, and 
falls to zero directly overhead. (The use of multiple RMCs 
with varying pointing directions will, however, smooth out 
sensitivity across the object scene.) 

Due to the fact that the RM incorporates multiple detectors, 
each with its own instrument response, the statistical signif- 
icance of an observation must include the variance of each 
detector's response: 




Pd(n) 2 - P d {n) 



SNR(n) « S{n) 



1 



P(n) 2 -P(n) 



(14) 



r\" \ OL) ' — ' l j / 

(15) 

where D is the total number of detectors. (The requirement 
that slat/slit width be equal to the detector diameter can be 
determined from this formula; unequal widths would decrease 
the profile contrast, inherently lowering the variance and thus 
the sensitivity.) Due to the addition of multiple response 
variances, the resulting sky sensitivity is more uniform, and 
due to the higher throughput, inherently better (approximately 
0.31) than the RMC. Because the PSF is distributed over a 
central peak and broad sidelobes (Fig. [3J, the RM sensitivity 
is, however, less than that of the coded aperture. Figure [4] 
shows the sensitivity curves for typical configurations of these 
three multiplexing instruments. 

VI. Reconstruction Technique 

The square-symmetric nature of P' n i n makes the system 
suitable for a wide variety of deconvolution techniques to 
solve for S, which have been developed across a broad 
spectrum of applications. For example, CLEAN |18| and 
the Maximum Entropy Method (MEM) lfl9l are widely used 
in radio interferometry, while the Algebraic Reconstruction 
Technique (ART) |20| and Maximum Likelihood Expectation 
Maximization (MLEM) |21| are commonly implemented in 
medical imaging. These traditional reconstruction techniques 
are unsuitable for the RM. CLEAN does not perform well 
when many sources are present due to the extended nature 
of the RM PSF, and it is unable to provide super-resolution. 
While we have previously shown imaging results using MEM 
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Fig. 4. Sensitivity plot for three multiplexing imagers, relative to Sy 



|[T5l . it is also unable to provide super-resolution or resolve 
multiple and complex sources. ART is subject to noise fluctu- 
ations due to its requirement that the reconstruction agrees ex- 
actly with the data. MLEM produces good results with super- 
resolution and suppression of noise, but is slow; the MLEM 
algorithm assumes that the data are Poisson distributed, and 
so deconvolution must take place directly from the temporal 
domain, where detector data can not be combined to speed up 
reconstruction. 

We have developed a new technique derived from Di- 
rect Demodulation (DDM) [22]. In simulations, DDM has 
been previously demonstrated to provide super-resolution with 
RMCs (23] and we have shown J24) the ability of DDM 
with an RM to achieve ~ 8x the geometric resolution and to 
resolve multiple and complex sources when little or no noise 
exists. For the reconstruction of actual measured data in the 
presence of background, we have modified the DDM routine 
and implemented a key step to compensate for noisy data. 
We refer to this technique as Noise-Compensating Algebraic 
Reconstruction (NCAR). 

The Gauss-Seidel iterative method is used to solve alge- 
braically for the object scene. The (k + 1) iteration of the 
reconstruction, f(n), is determined by E5ll 



f ik+1) (n) 



P 



- ]T P'{m,n)f k {m)) , (16) 



with f°(n) = 0. Gauss-Seidel uses the results of calculations 
from the same iteration, and so to prevent pixel bias, it 
is prudent to randomize the order of m every iteration. A 
positivity condition is enforced as a physical constraint on 
each pixel of f^ k+1 \n). 

A normalized parameter f3 may be defined, which sums 
over all pixels in the residual to examine the agreement of 



ft 0.5 : 




(a) 




(b) 



Fig. 5. Simulation of a single point source centered in the FOV, showing the 
reconstructed image (with actual source location designated by a + symbol) 
and corresponding normalized residual summation, 0, plotted against iteration 
number using (a) Direct Demo dula tion, and (b) with the noise compensation 
additive, R[cr(n)], given in Eq. |18| In (b), the vertical dashed line corresponds 
to ft, the iteration at which the maximum agreement has been reached, and 
the horizontal dashed line corresponds to the residual convergence value given 
the noisy data. Note the contrast to (a), which converges indefinitely to zero 
(i.e. perfect agreement). 
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the reconstruction with the cross-correlation, 



P 



(fc+i) 



1 



2J 



J2P'(n',n)f( k+1 \n)-C(n') 



(17) 



In Fig. [5^, this residual summation is plotted against the 
iteration number for a sample reconstruction. The successive 
reconstructions converge to the data (i.e. /3 shrinks) for an 
indefinite number of iterations. This agreement, however, 
includes any fluctuations from noise within the data, and so 
this method causes spurious peaks and poor location accuracy 
of the true source, as seen in the reconstruction, when the 
signal-to-noise ratio (SNR) is low. 

This "noise amplification" is the primary problem in image 
reconstruction using an algebraic technique. Regularization 
methods (e.g. smoothing, pixel-to-pixel variation constraints) 
may be employed, but this regularized image is still forced to 
agreement with the data, C(n). The NCAR technique allows 
for deviation from the noisy data, by replacing C(n) in Eq. 
16] 



with C( fc+1 )(n), such that 
C {k+1 \n) = C{n) 



i?( fc+1 )[ f j(n)]. 



(18) 



The function R[a] provides a random number from a normal 
distribution with average zero and standard deviation a, where 
<r(n) is the error in each pixel of the cross-correlation. If the 
error is assumed to be only from Poisson uncertainty, it can 
be shown that 



<r(n) = k jy ] 0(m)P(m, n) 2 . 



(19) 



The reconstruction and residual summation, f3, are plotted 
for NCAR in Fig. [5j?. The agreement between reconstruc- 
tion and data fluctuates about an envelope showing initial 
convergence followed by a leveling off beginning at some 
iteration k. Once this minimum f3 is reached, the image has 
converged as well as possible to the data given the noise. 
Each successive iteration still contains spurious peaks arising 
from noise fluctuations, but because the component R[a(n)] 
has been added to randomize the noise, the amplitude and 
location of these peaks are unique to each iteration, while the 
iterative images share the true sources. By now averaging over 
all the images after iteration k, the spurious random peaks are 
suppressed while the true source distributions remain: 

f(k+i)(n) + (k - K)F k (n) 



F^ k+1 \n) 



k+1 > k. (20) 



k + 1 - k 

F^ k+1 ^(n) is the final reconstruction shown in Fig. [s^. Using 
NCAR, statistical uncertainty in the data is manifest as an 
increase in the size of the reconstructed source in the form of 
fuzzy sidelobes. In cases where the actual size of the source is 
desired, NCAR may not be desirable; but with large noise in 
the data, an algebraic technique alone will produce spurious 
peaks and may not represent the true source size anyway. 
When it is acceptable that source size in the reconstructed 
image relates to uncertainty in the data, however, NCAR 
produces image reconstructions with a significant reduction in 
noise fluctuations compared to results using MEM, CLEAN, 
or algebraic solutions, as will be demonstrated in detail in Sec. 
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Fig. 6. Contour plot of the relative sensitivity of LaBRAT's FOV, with 
the parameters used for the simulations described in Sections |VII| and |VIII| 
Values are normalized to S{T/B) 1 / 2 . Axes are in decrees. 



VII. Simulation 

Monte Carlo simulations were performed to test and verify 
the NCAR technique. The software was custom written for this 
project in IDL. The instrument response and cross-correlation 
calculated 



images are calculated as described in Sec. IV Since the 
instrument response assumes far-field imaging, sources are 
simulated at an "infinite" distance from the detector; i.e. pho- 
tons from the same source arrive at equal angles of incidence. 

The geometry of the simulated RM reproduces that of a 
laboratory prototype, the Lanthanum Bromide-based Rotat- 
ing Aperture Telescope (LaBRAT), developed at LSU |15|. 
LaBRAT's detection plane is composed of nineteen cylindrical 
Cerium-doped Lanthanum Bromide (Ce:LaBr3) scintillators, 
each 3.8 cm diameter x 2.5 cm thick. They are positioned 
in a hexagonal layout and coupled to 3.8 cm PMTs. The 
mask is a grid of eight 3.8 cm x 61 cm x 1.9 cm thick 
lead slats spaced 3.8 cm apart and sitting ~1.2m away from 
the detectors. This arrangement defines a 13.8° diameter FOV 
and a 1.9° FWHM geometric angular resolution. As a result 
of the 19 detectors used in LaBRAT, the sensitivity is fairly 
uniform. As seen in Fig. [6] where white indicates regions of 
the object scene that are unmodulated from the viewpoint of 
one of the detectors (i.e. has less imaging sensitivity), the 
overall relative sensitivity is close to 0.31 over the entire scene. 
The pattern in the figure thus mimics the detector layout. 
Additional dithering of the telescope orientation about the 
vertical (as would occur in a balloon-borne instrument) would 
result in additional uniformity across the FOV. 



In Sec. VIII the NCAR technique is compared to its precur- 



sor, DDM, to demonstrate the effect of noise compensation. 
Additionally, we include a comparison to MEM and CLEAN, 
because of their wide use in astrophysical applications. For 
these simulations, the FOV is divided into 12x12 arcmin bins, 
totaling 4900 pixels for a single image. Each of the results is 
that of a 30-minute simulated exposure. NCAR, DDM, and 
CLEAN are run for 10k iterations, and MEM for 50k - 100k 
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Fig. 7. Monte Carlo Simulation results for a single point source and no 
background. Object scene is shown at upper left, raw cross-correlation image 
(Eq. [12) is at upper right, and results of MEM, CLEAN, DDM and NCAR 
are shown below. All reconstruction techniques perform reasonably well. 
(Axes in degrees) 



Fig. 8. Monte Carlo Simulation results for a single point source with 
background. Note that, in the presence of background, the MEM, CLEAN, 
and DDM reconstructions contain spurious point-like sources and the image 
does not convey the uncertainty of the location of the true source. NCAR, 
however, smooths over the noise contributions, and places an uncertainty on 
the true source that is related to the SNR of the measurement. (Axes in 
degrees) 



to assure convergence. 

VIII. Results 

We first observed how a single point source with a measured 
rate of 10 Hz (over all detectors) is imaged in the presence and 
absence of background. When no background exists (Fig. 17), 
all deconvolution techniques reconstruct the scene reasonably 
well, although spurious point sources are faintly visible in the 
MEM reconstruction. 

When a background rate of 500 Hz is included (Fig. [8), 
the MEM, CLEAN, and DDM reconstructions exhibit noise 
fluctuations that are manifest as spurious peaks in the image. 
The true source itself, however, remains relatively unaffected. 
The NCAR image shows the presence of diffuse background 
features, and the true source peak has broadened. Due to the 



nature of NCAR, the width of a point source reconstruction 
narrows as the SNR grows. This can be a desirable feature of 
an image, where a "blur" is typically considered an uncertainty 
on source location. As we continue to show, this blurring 
characteristic does a good job of removing spurious noise 
fluctuations. 

In Figs. [9p0] object scenes with two sources, each with a 
measured rate of 50 Hz, and a background contribution of 100 
Hz, are simulated. The background is kept relatively low here 
to demonstrate the intrinsic ability of NCAR to achieve super- 
resolution. In Fig. [9] the sources are separated by 3°, greater 
than the 1.9° geometric resolution of the instrument. NCAR 
produces the result that is the most free of spurious peaks and 
other degrading artifacts. Due to the broadened nature of the 
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Fig. 9. Monte Carlo Simulation results for two equal intensity sources 
separated by 3°. Statistical and algebraic reconstruction techniques are both 
capable of resolving two sources that are at an angular separation greater than 
the geometric resolution defined by the instrument. (Axes in degrees) 



Fig. 10. Monte Carlo Simulation results for two equal intensity sources 
separated by 1°. At angular separations less than the geometric resolution 
of the instrument, MEM reconstructs a single elongated source and CLEAN 
sees only a single source, while DDM and NCAR are able to fully resolve 
the two sources. (Axes in degrees) 



RM PSF, CLEAN in particular shows difficulty in resolving 
sources at this separation. In Fig. 10 the sources are moved 
to a separation of 1°, about half the intrinsic resolution of the 
instrument. MEM reconstructs an elongated structure, while 
CLEAN reconstructs only a single point source. Additional 
simulations have confirmed that MEM and CLEAN are both 
limited by the geometric resolution of the instrument. DDM 
and NCAR resolve both sources and provide the desired super- 
resolution. This result has also been previously demonstrated 
with measured laboratory data, resulting in a resolving power 
of 20', or ~ 8x the geometric resolution of the instrument 
[26 1 . Again, NCAR produces the most accurate reconstruction 
of the object scene. 

NCAR also provides image reconstructions for more com- 
plicated source distributions. Fig. 11 shows the results for 
5 sources, each with a measured source rate of 50 Hz and 
equivalent background rate. MEM is capable of resolving 



some, but not all of the sources entirely, while CLEAN 
performs slightly better. DDM suffers from spurious peaks 
near the true sources. NCAR, however, reconstructs each of 
the sources well, with no visible noise, but instead, a slight 
blurring. 



A line source distribution (Fig. [12) is also simulated. To 
examine the ability of the reconstruction techniques to resolve 
an extended source, the background rate is set to zero, and 
the line source has a total measured rate of 4.5 kHz. MEM 
performs poorly with a significant spurious reconstruction to 
the left due to the interference of sidelobes. The CLEAN 
result has no significant spurious reconstructions, but the line 
is poorly resolved. DDM suffers from some spurious peaks 
arising from noise and interfering PSFs from the complex 
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Fig. 11. Monte Carlo Simulation results for five equal intensity sources. Fig. 12. Monte Carlo Simulation results for a line source distribution. 

Even with little background in the measurement, both MEM and DDM suffer The PSF interference patterns cause MEM to perform poorly with a complex 

from noise artifacts in the reconstruction due to the more complex object source, while DDM continues to suffer from from noise artifacts. NCAR 

scene. NCAR, however, compensates for the noise and produces an image removes most of the spurious sources. (Axes in degrees) 
better representative of the object scene. (Axes in degrees) 



nature of the object scenes. NCAR provides a reconstruction 
of the line similar to DDM with most of the spurious peaks 
smoothed out. 

The "broken" lines observed in the DDM and NCAR 
reconstructions in Fig. [T2] are not due to sensitivity variations. 
As seen in Fig. [6] the sensitivity is relatively uniform, and 
the low-sensitivity sky regions are uncorrelated with the line 
breaks. Rather, this is a consequence of the statistical noise 
in the image causing the deconvolution to converge on an 
imperfect solution to the data. This effect is seen to increase for 
more complex object scenes (i.e. the number of non-trivial sky 
bin values increases), and is an inherent property of the RM 
response, not the deconvolution algorithm. The effect could 
be reduced by increasing the exposure time (thus improving 
SNR). 

In order to observe two weak sources imaged simultane- 



ously, a 25 and 15 Hz source, in the presence of 500 Hz 
background, are simulated in Fig. [13] MEM performs poorly, 
showing spurious reconstructions indistinguishable from the 
weaker of the two sources. CLEAN and DDM perform the 
best, with two sources clearly reconstructed with relative fluxes 
apparent. The "blurring" characteristic of NCAR, however, 
means that the weaker of the two sources is blurred to a greater 
degree, suppressing its apparent brightness further. For higher 
SNR, this discrepancy disappears. 

IX. Discussion 

For imaging modalities with 50% transmission, the coded 
aperture offers the maximum sensitivity possible. By moving 
to the temporal domain for photon modulation, the RM gains a 
unique advantage in that the measurement vector is effectively 
a continuous function, and so finer sub-sampling allows one 
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roughly S9 > a/apL ps 2.6°. In other words, at hard X- 
ray energies and above, a coded aperture has a sensitivity 
advantage over an RM but cannot take advantage of its full 
potential for excellent angular resolution. At high energies, 
the super-resolution capability of an RM with NCAR allows 
it to provide good angular resolution although there is a 
penalty of 1.6r/ 2 in the sensitivity compared to the coded 
aperture. At lower energies, the RM can utilize its super- 
resolution capability, while the coded aperture is restricted to 
the geometric angular resolution. 

When compared to the RMC, the RM maintains a sensi- 
tivity advantage of 0.31/0.15 w 2, due to its higher mask 
transmission. Additionally, the moderate spatial resolution of 
the RM detection plane increases the number of independent 
measurements of the object scene, contributing to enhanced fi- 
delity of the reconstructed image and smoothing the sensitivity 
curve across the sky, and the absence of a thick secondary grid 
reduces instrument weight. 

An algebraic technique such as NCAR is necessary to 
exploit the super-resolution capability of the RM while si- 
multaneously suppressing noise fluctuations. NCAR will not, 
however, be suitable for all scenarios. As shown, its weakest 
performance is in resolving two weak sources with different 
measured rates. Because of the inherently low SNR, the peaks 
are "blurred" to different degrees, suppressing the weaker 
of the two and thus not accurately depicting relative source 
strength. In the other cases presented, however, NCAR per- 
forms quite well, reconstructing the object scene and suppress- 
ing noise fluctuations that typically plague algebraic solutions. 
In scenarios where the size of the source is of importance, 
NCAR will not produce good results with low SNR. If the 
ultimate desire is to survey and locate true sources, however, 
NCAR generates images free of spurious reconstructions with 
a visual representation of the locational uncertainty of a 
particular source measurement. 



Fig. 13. Monte Carlo Simulation results for weak source (lower right) in 
the presence of a stronger source (upper left). (Axes in degrees) 



to achieve super-resolution with an appropriate reconstruction 
algorithm. The RM does, however, suffer in sensitivity due to 
its non-ideal response. 

At the geometric angular resolution 69 = a/L (i.e., 
with super-resolution factor r] = 1), a coded aperture has 
an intrinsic sensitivity advantage over an RM by a factor 
0.50/0.31 « 1.6 (Eqs. [P] [BJ. The width/diameter b of the 
detector pixel is a fraction a of the mask element size a: 



b = a/a (typically, a = 1 for an RM and 2-4 for a coded 
aperture). For good efficiency at high energies, the detector 
width and thickness should be at least equal to the photon 
interaction length, 1/crp, implying that 



ap 



> 



LS8 



(21) 



i.e., there is a minimum ap (or equivalently a maximum 
energy) at which a particular angular resolution is achievable. 
For example, at 511 keV, with L = 1 m and a = 2, the 
angular resolution of a coded aperture imager is limited to 



X. Conclusion 

A rotating modulator is an instrument capable of imaging 
at hard x-ray and gamma-ray energies. The instrument is 
composed of a single grid of slats above a collection of circular 
non-imaging detector elements. While the instrument response 
for the RM is non-ideal (the sensitivity is approximately 62% 
of that of an equivalent coded aperture), the near-continuous 
nature of the observed count profiles enables an analysis 
which can go beyond the instrinsic resolution as defined 
by the geometry, whereas the coded aperture cannot. When 
compared to the RMC, the RM is able to achieve higher 
and more uniform sensitivity, does not suffer from mechanical 
collimation at high energies, and has lower overall weight. 

We have found that algebraic solutions are the only recon- 
struction techniques that can both achieve super-resolution and 
perform relatively fast. Algebraic methods, however, have in 
the past suffered from poor reconstructions due to the require- 
ment of strict agreement with the noisy data. We have devel- 
oped a novel technique, NCAR, based on an algebraic solution 
and non-linear physical constraints. A noise-compensation 
factor derived directly from the Poisson uncertainty in the data 
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is added to a Gauss-Seidel iteration. By performing a running 
average on the successive results after a specific agreement 
criterion has been met, the spurious peaks arising from noise 
in the data are effectively suppressed, while the true sources 
remain. The uncertainty of the true source location due to noise 
in the data is manifest as a blurring of the scene rather than 
spurious peaks elsewhere in the image. 

When compared to other common deconvolution algorithms 
(MEM, CLEAN, and DDM), NCAR provides higher fidelity 
in most cases. Furthermore, the technique is general enough 
that it could be used in imaging and other applications where 
deconvolution is required in the presence of large background. 
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